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Abstract 

Almost  no  data  exists  in  the  literature  on  thermochemical  properties-  enthalpy  of  formation, 
entropy  or  heat  capacities  -  for  chlorinated  sulfur  hydrocarbons  and  oxygenated  intermediates  in 
atmospheric  or  thermal  oxidation  reactions.  Even  less  kinetic  data  exists  on  the  elementary 
processes  involved  in  these  reaction  systems.  This  study  describes  a  first  attempt  to  generate  a 
thermochemical  property  data  set  and  an  elementary  reaction  mechanism  for  pyrolysis  and 
oxidation  of  diethyl  sulfide  and  dichloro-diethyl  sulfide. 

Standard  enthalpy  of  formation, A//°29g,  entropy,  S°9 8,  and  heat  capacities  CAT),  are 

determined  for  stable  molecules,  intermediates,  and  transition  states  in  oxidation  of  sulfur 
containing  hydrocarbons  using  DFT  (Density  Functional)  calculations  with  the  B3LYP  and  the 
KM-LYP  functionals.  Ab  initio  calculations,  using  the  CBS-Q  composite  method  are  also  used. 
Benson-type  groups  for  group  additivity  estimation  of  oxygenated  sulfur  hydrocarbons  are 
constructed. 

Enthalpies  of  formation  (A/jYVa)  are  determined  using  the  Density  Functional  and  CBSQ 
methods  and  working  reaction  analysis  with  isodesmic  or  group  balance  constraints.  Entropy  (S) 
and  heat  capacity  ( Cp(7 )),  values  from  vibrational,  translational,  and  external  rotational 
contributions  are  calculated  using  the  rigid-rotor-harmonic-oscillator  approximation,  based  on 
the  vibration  frequencies  and  structures  obtained  from  the  density  functional  studies.  The 
contribution  to  S  and  Cp  (T)  from  analysis  on  the  internal  rotors  is  incorporated  into  the  analysis 
instead  of  the  more  conventional  use  of  torsion  frequencies.  This  is  an  important  consideration  in 
the  thermochemical  and  kinetic  analysis;  because  of  the  loss  of  internal  rotors  in  transition  states. 
The  important,  (low  energy  barrier)  paths  for  unimolecular  dissociation  of  the  parent  sulfur 
chlorocarbons  are  HC1  elimination,  retro-ene,  and  carbon  sulfur  bond  cleavage  reactions.  These 
are  however,  all  above  50  kcal/mol  in  barrier. 

Details  on  the  thermochemistry  and  kinetic  parameters  for  thermal  reactions  are  presented 
herein.  A  detailed  pyrolysis  and  oxidation  mechanism  is  constructed  using  hydrocarbon  and 
chlorocarbon  oxidation  using  nitrogen-oxygen  mechanisms  from  the  literature,  and  sulfur- 
hydrocarbon  reactions  from  this  study.  A  detailed  chemical  kinetic  mechanism  for  HD  pyrolysis 
and  oxidation  has  been  created  using  kinetics  and  thennodynamic  properties  generated  in  this 
work  combined  with  published  sulfur,  chlorine,  nitrogen,  and  hydrocarbon  kinetics.  Sensitivity 
analysis  is  used  to  detennine  important  HD  destruction  pathways. 
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Introduction 


Bis(2-chloroethyl)  sulfide,  (ClCfifiCfifi^S  is  a  main  component  of  the  chemical  warfare  agent 
HD  or  mustard  gas,  as  has  been  recognized  since  the  late  1 800's.  It  has  been  manufactured  in 
quantity  and  relatively  pure  form  since  the  early  191 0's.  More  recently,  within  the  last  decade  or 
so,  chemical  warfare  agents  are  being  destroyed  due  to  their  potentially  harmful  effects  to 
civilian  populations.  The  HD  (bis(2-chloroethyl)  sulfide)  is  a  liquid,  that  has  a  relatively  low 
vapor  pressure,  relative  to  phosgene  (Cl2C=0  or  hydrogen  cyanide  (HCN).  The  low  volatility  is 
due  to  its'  moderately  high  molecular  weight,  it  boils  near  220  C.  The  high  toxicity  of  HD, 
similar  to  that  of  phosgene  and  HCN,  make  experimental  studies  on  this  species  difficult.  In 
effect,  there  are  very  limited  studies  on  it  and  other  agents  ’  ,  and  their  reliability  has  been 
critisized3. 

Mustard  gas  was  manufactured  in  large  quantities  during  World  Wars  I  and  II,  but  has  not  been 
manufactured  on  an  industrial  basis  in  the  United  States  since  1968.  Stockpiles  of  mustard  gas 
are  stored  at  seven  locations  scattered  around  the  United  States.  Incineration  has  been  selected 
by  the  United  States  as  one  of  the  methods  for  destroying  highly  toxic  chemical  agents  and 
munitions  contained  within  the  chemical  weapons  stockpile.  As  part  of  the  Chemical  Stockpile 
Disposal  Program  mandated  by  Congress,  the  U.S.  Army  currently  uses  the  baseline  incineration 
system  for  destruction  of  munitions  and  bulk  agents  containing  mustard  gas.  Incineration 
technology  has  been  used  at  two  sites,  and  is  under  construction/systemization  at  three  additional 
sites.  The  scientific  and  engineering  aspects  of  chemical  weapons  incineration  were  reviewed  by 
Gouldin  &  Fisher4.  A  simulation  took  for  analyzing  chemical  demilitarization  incinerators  is 
being  developed  to  help  reduce  the  time,  cost  and  technical  risk  of  destroying  the  stockpile 
(Denison).  Such  a  tool  requires  a  detailed  model  of  the  chemical  processes  occurring  within  the 
incineration  equipment.  A  detailed  chemical  kinetic  mechanism  for  combustion  of  the  nerve 
agent  GB  has  been  developed6,  but  no  kinetics  or  thermochemistry  existed  previously  for 
modeling  combustion  of  HD  or  the  nerve  agent  VX. 

Combustion  or  thermal  reaction  of  HD  in  the  presence  of  oxygen  to  convert  it  into  the  products 
CO2,  HC1,  SO2  (or  H2SO4)  and  H20  is  one  method  that  can  be  used  effectively  to  destroy  the 
mustard  agent.  The  combustion  or  oxidation  process,  however  needs  to  completely  destroy  the 
agent  and  its  immediate  decomposition  products.  The  thermal  reactions  of  mustard  or  any 
hydrocarbon  sulfide,  include  pyrolysis  (unimolecular  decomposition),  as  well  as  reactions  with 
the  radical  pool  to  effect  hydrocarbon  and  sulfur  radicals  of  the  bis(2-chloroethyl)  sulfide.  These 
radicals  then  react  with  oxygen  molecules  and  oxygen-containing  radical  species  (O,  OH,  H02, 
CH3O,  etc.),  eventually  converting  to  lower  energy  products,  as  noted  above. 

Sulfur  species  are  also  present  in  many  fuels,  including  coals,  heating  oil,  aircraft  fuels,  and 
diesel  and  gasoline.  In  all  these  cases  the  source  of  fuel  sulfur  is  contact  of  the  elemental  forms 
of  sulfur  with  the  natural  underground  deposits,  where  sulfur  becomes  incorporated  into  the  gas 
and  crude  oil  stocks.  A  significant  fraction  of  the  sulfur  in  these  liquid  fuels  is  removed  prior  to 
use  (combustion)  in  order  to  prevent  pollution  of  the  air,  formation  of  S02  and  subsequent 
formation  of  particulates  and  acid  rain.  It  is  more  difficult  to  remove  the  sulfur  from  coals;  but 
this  is  performed  in  many  cases  as  well.  Sulfur  does  persist  in  these  cleaned  fuels  and  oils  that 
are  treated  for  its  removal  to  some  extent  and  the  sulfur  species  are  subsequently  combusted  in 
their  use. 
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In  plants  and  animals,  sulfur  occurs  in  various  proteins,  and  it  is  one  of  the  10  most  abundant 
elements  in  the  human  body.7  Sulfur  compounds  are  present  and  of  interest  in  biological 
systems,  in  atmospheric  chemistry,  and  in  environmental  science.  The  known  involvement  of 
sulfur-centered  radicals  in  these  biological  systems  suggests  that  knowledge  on  the  fundamental 
thermodynamic  properties  of  the  sulfur  moieties  would  be  valuable  in  understanding  its  effects 
on  these  systems.8 

The  knowledge  and  understanding  the  oxidation  and  thermal  reactivity  of  these  sulfur  species 
requires  their  thennodynamic  properties,  bond  energies  and  kinetic  parameters  for  unimolecular 
reactions  to  be  known. 

The  use  of  Density  Functional  Theory,  B3LYP/6-31g(d,p),  with  isodesmic  working  reactions  for 
enthalpy  of  formation  of  sulfur  hydrocarbons  is  tested  in  this  study  using  a  set  of  known  sulfur 
hydrocarbons.  Satisfactory  results  (below)  indicate  that  we  can  use  this  computational  chemistry 
method  to  determine  thennochemical  properties  and  kinetics.  This  part  I  of  our  study  focuses  on 
unimolecular  (thermal)  reactions  of  the  mustard  and  on  its  thermochemistry  for  evaluation  of 
initial  abstraction  reactions  that  result  in  the  initial  radical  products.  These  radical  intermediates 
react,  as  noted  above,  with  oxygen  species  and  by  further  unimolecular  decomposition  processes. 

Thermodynamic  properties  for  reactants,  transition  states,  and  products  from  unimolecular 
dissociations  of  CH3CH2SCH2CH3,  CH3CH2SCH2CH2CI,  and  CFFClCFFSCFFCFbCl  and  for 
their  radicals  resulting  from  bond  cleavage  reactions  are  analyzed.  Af//°29s  for  the  molecules  and 
radicals  are  determined  using  isodesmic  reaction  analysis  at  the  B3LYP/6-31G(d,p)  level,  with 
S° 298  and  Cp(T)  determined  using  geometric  parameters  and  vibrational  frequencies  obtained  at 
this  same  level  of  theory.  Potential  barriers  for  the  internal  rotor  potentials  are  also  calculated  at 
the  B3LYP/6-31G(d,p)  level,  and  the  hindered  rotation  contributions  to  S°29 8  and  Cp(T)  are 
calculated.  A  more  limited  number  of  studies,  specifically  for  transition  states  are  also  performed 
at  KMLYP/6-31  lG(d,p),  CBS-Q,  and  CBS-QB3  levels  of  theory,  with  the  CBS-QB3  barriers 
recommended.  Hydrogen  bond  dissociation  energies  of  CH3SH,  CH3CH2SH,  CH3SCH3, 
CH3CH2SCH2CH3,  and  CH3CH2SCH2CH2CI  are  calculated  and  compared  with  literature  values 
where  available.  Four  center  (Retro-ene  like)  reactions  and  C — S  bond  dissociations  are 
important  initial  decomposition  channels  in  the  hydrocarbon-sulfur  moieties. 

Calculation  Method 

Geometries  and  Vibration  Frequencies  The  initial  structures  of  reactants,  transition  states  and 
products  are  determined  using  ROHF  or  UHF/PM3  in  MOPAC,9  followed  by  optimization  and 
vibrational  frequency  calculation  at  B3LYP/6-31G(d,p)  level  using  GAUSSIAN  98. 10  Transition 
state  geometries  are  identified  by  the  existence  of  only  one  imaginary  frequency,  structure 
information,  and  the  transition  state  theory  reaction  coordinate  vibration  information.  Zero-point 
vibrational  energies  are  scaled  by  0.9806  as  recommended  by  Scott  and  Radom.11  Geometries  of 
transition  states  optimized  at  B3LYP/6-31G(d,p)  are  shown  in  Figure  1. 

Enthalpies  of  Formation  We  evaluated  enthalpies  of  formation  (Af/ /°29x)  at  B3LYP/6-31G(d,p) 
level  on  several  sulfur  hydrocarbons  and  corresponding  radicals  with  a  method  of  isodesmic 
reaction  and  show  that  the  B3LYP/6-31G(d,p)  values  are  in  agreement  with  accepted  literature 
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values  in  Table  l.The  evaluated  At-fl°29H  for  the  reference  molecules  and  radicals  in  the  isodesmic 
reactions  are  given  in  Table  2. 

Enthalpies  of  formation  for  reactants  and  products  are  calculated  at  B3LYP/6-31G(d,p)  level. 
The  method  of  isodesmic  reactions  for  OEECTESH,  CH3OHSH,  CH3CH2SCH2OH2, 
CH3CH2SOHCH3,  CH3CH2SCH2CH2C1,  CH3CH2SCH2C«HC1,  and  CH3CH2SC«HCH2C1  is 
used  to  determine  for  accurate  enthalpies  of  formation  at  B3LYP/6-31G(d,p)  and  recently 
published  KMLYP  density  functional  method  with  the  /6-311G(d,p)  basis  set.12  The  evaluated 
reaction  enthalpies  and  enthalpies  of  formation  in  the  isodesmic  reactions  are  listed  in  Table  3. 


C*H2CH2SH  +  CH3SH  <=>  CH3CH2SH  +  C*H2SH . (IR1) 

CH3C*HSH  +  CH3SH  <=>  CH3CH2SH  +  C*H2SH . (IR2) 

CH3CH2SCH2C*H2  +  CH3CH2SH  CH3CH2SCH2CH3  +  C*H2CH2SH . (IR3) 

CH3CH2SC*HCH3  +  CH3CH0SH  <=>  CH3CH2SCH2CH3  +  CH3C*HSH  . (IR4) 

CH3CH2SCH2CH2C1  +  CH3CH3  CH3CH2SCH2CH3  +  CH3CH2C1  . (IR5) 


CH3CH2SCH2C*HC1  +  CH3CH2OH  <=>  CH3CH2SCH2CH2C1  +  C*H2CH2OH  (IR6a) 

CH3CH2SCH2C*HC1  +  CH3OCH3  CH3CH2SCH2CH2C1  +  CH3OC*H2  ,.(IR6b) 

CH3CH2SC*HCH2C1  +  CH3CH2OH  <=>  CH3CH2SCH2CH2C1  +  CH3C*HOH  (IR7a) 

CH3CH2SC*HCH2C1  +  CH3OCH3  «■  CH3CH2SCH2CH2C1  +  CH3OC*H2  ,.(IR7b) 

Enthalpies  of  formation  for  transition  states  Although  B3LYP  has  achieved  good  success  in 
predicting  thermochemistry  of  molecules,  several  studies  report  that  B3LYP  is  less  accurate  for 
predicting  barriers  and  has  a  tendency  to  underestimate  barriers.12'15  We  performed  further 
calculations  for  barriers  with  four  different  methods  (B3LYP/6-31G(d,p),  KMLYP/6- 
311G(d,p)12,  CBS-Q16'18,  and  CBS-QB319,20).  Comparison  of  barriers  is  given  in  Table  4.  The 
barrier  is  calculated  relative  to  the  reactants.  The  complete  basis  set  (CBS-QB3)  method  of 
Montgomery  et  al. 19,20  for  computing  accurate  energies  is  used  to  detennine  the  activation 
energies  used  in  our  kinetic  analysis. 

Entropy  and  Heat  Capacity  The  contributions  of  external  rotations,  translations,  and  vibrations 
to  entropies  and  heat  capacities  are  calculated  from  scaled  vibration  frequencies  and  moments  of 
inertia  for  the  optimized  B3LYP/6-31G(d,p)  structures.  The  number  of  optical  isomers  and  spin 
degeneracy  of  unpaired  electrons  are  also  incorporated.  Contributions  from  hindered  internal 
rotation  for  S  and  Cp(  T)  are  determined  using  direct  integration  over  energy  levels  of  the 
intramolecular  rotational  potential  curves.  A  program,  ROTATOR  ,  is  used  for  calculation  of 
the  energy  levels.  This  technique  employs  expansion  of  the  hindrance  potential  in  the  Fourier 
series  (El),  calculation  of  the  Hamiltonian  matrix  on  the  basis  of  wave  functions  of  the  free 
internal  rotor,  and  subsequent  calculation  of  energy  levels  by  direct  diagonalization  of  the 
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Hamiltonian  matrix.-  The  torsional  potential  calculated  at  discrete  torsion  angles  is 
represented  by  a  truncated  Fourier  series: 

V(&)  =  oo  +  2  a,- cos (i0)  +  X  bt sin(/0)  i=  1,2,3,  ...  (El) 

Values  of  the  coefficients  (ao,  o,  and  b,)  are  calculated  to  provide  the  minimum  and  maximum  of 
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the  torsional  potentials  with  allowance  of  a  shift  of  the  theoretical  extreme  angular  positions-  ' 
The  energy  levels  are  used  to  determine  partition  coefficients  and  their  contributions  to  S  and 
Cp{T)  through  relationships  from  statistical  mechanics. 

High-Pressure  Limit  Pre-exponential  Factors  and  Rate  Constants  For  the  reactions  where 
thennochemical  properties  of  transition  states  are  calculated  by  ab  initio  or  density  functional 
methods,  Ws  are  fit  by  three  parameters  A,  n,  and  Ea  over  temperature  range  from  298  to 
2000K,  k^,  =  A(T)n  exp  {-Ea  /RT).  Entropy  differences  between  reactant  and  transition  state  are 
used  to  determine  the  pre-exponential  factor,  A,  via  canonical  Transition  State  Theory  (TST) 

A  =  (khT/hp)cxp(ASf/R),  Ea  =  A Hx 

Where  hp  is  the  Planck  constant  and  kb  is  the  Boltzmann  constant.  Treatment  of  the  internal 
rotors  for  S  and  Cp(  T)  is  important  here  because  these  internal  rotors  are  often  lost  in  the  cyclic 
transition  state  structures. 

Results  and  Discussion 

Rotational  Barriers  Potential  barriers  for  internal  rotations  of  CH3CH2SCH2CH2CI  and 
CH2CICH2SCH2CH2CI  are  calculated  in  the  B3LYP/6-31G(d,p)  level  and  shown  in  Figures  2 
and  3,  respectively.  The  potential  energies  are  calculated  as  a  function  of  the  dihedral  angle  by 
varying  the  torsion  angle  in  30°  intervals  and  allowing  other  parameters  to  be  optimized.  The 
barriers  for  internal  rotations  are  calculated  from  the  differences  between  the  total  energy  of  each 
conformation  and  that  of  the  most  stable  conformer  at  0  K,  where  the  zero  point  vibrational 
energy  (ZPVE)  and  thermal  correction  to  298  K  are  not  included.  Total  energies  at  0  K  and 
calculated  internal  rotation  barriers  versus  the  dihedral  angle  are  presented  in  the  coefficients  of 
the  Fourier  expansion  components,  a,,  and  b,  in  equation  El. 

Entropy  (5°?qs)  and  Heat  Capacity  ( Cp(T ))  S°2 98  and  Cp(  7)’s  are  calculated  based  on  scaled 
vibration  frequencies  and  moments  of  inertia  of  the  optimized  B3LYP/6-31G(d,p)  structures. 
The  calculation  results  are  summarized  in  Table  5.  TVR  represents  the  sum  of  the  contributions 
from  translation,  vibrations  and  external  rotations  for  5°(298)  and  Cp(T)A.  Symmetry,  number  of 
optical  isomers  and  electronic  spin  are  incorporated  in  estimation  of  S° (298)  as  described  in  Table 
5.  Torsion  frequencies  are  omitted  in  these  calculations,  instead,  contributions  from  internal 
rotation  for  5°(298)  and  Cp(T)' s  are  determined  using  direct  integration  over  energy  levels  of  the 
intramolecular  rotational  potential  curves14'16  and  noted  in  Table  5. 

Hydrogen  bond  dissociation  energies  in  CH3SH,  CH3CH?SH,  CH3SCH3,  CH3CH?SCH?CH3,  and 
CH3CH2SCH2CH2CI  Hydrogen  bond  dissociation  energies  (BDEs)  for  the  C — H  and  S — H 
bonds  in  CH3SH,  CH3CH2SH,  CH3SCH3,  CH3CH2SCH2CH3,  and  CH3CH2SCH2CH2C1  are 
presented  in  Table  6.  Our  calculation  results  at  B3LYP/6-31G(d,p)  level  of  theory  are  compared 
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with  literature  values  when  available.  If  not,  the  values  of  KMLYP  density  functional  method 
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with  the  /6-31  lG(d,p)  basis  set  are  compared.  For  alkyl  RSH  species  Lias  et  al7  and  Griller  et 
al.  agreed  that  BDE(S — H)  was  effectively  independent  of  the  structure  of  the  R  group.  This 
was  based  on  values  of  BDE  derived  from  the  electron  affinities  of  the  RS  radicals  and  the 
A//add°  values  of  the  molecules.  They  recommended  a  value  of  87.00  ±  2.15  kcal/mole  '  and 
88.43  ±  1.91  kcal/mol26,  respectively,  for  BDE(S — H).  Amstrong6  also  reported  BDE(S — H) 
value  for  all  the  RSH  species  is  87.47  ±  0.96  kcal/mol.  Our  calculation  values  are  ~  89  kcal/mol 
for  CH3S — H  and  CH3CH2S — H.  We  determine  the  BDEs  of  H— CH2SH  (96.3)  and  H— 
CH2SCH3  (95.4)  are  about  2  kcal/mol  higher  than  literature  values.  (93.9  and  93.3  kcal/mol, 
respectively.25,27'29)  The  corresponding  hydrogen  BDEs  for  CH3CH2SH,  CH3CH2SCH2CH3,  and 
CH3CH2SCH2CH2C1  are  reported  in  Table  6. 

Unimolecular  Dissociation  Reactions  of  CH3CH2SCH2CH3,  CH3CH?SCH2CH2C1,  and 
CH2C1CH2SCH2CH2C1  Potential  energy  diagrams  for  unimolecular  dissociations  of 
CH3CH2SCH2CH3,  CH3CH2SCH2CH2C1,  and  CH2C1CH2SCH2CH2C1  are  illustrated  in  Figure  4. 

Retro-ene  Reaction  CH3CH2SCH2CH3  can  dissociate  to  CH3CH2SH  +  C2H4  via  a  four-center 
hydrogen  shift  transition  state  (TS1).  This  reaction,  when  involving  a  six-member  ring,  is  often 
tenned  a  retro-ene  reaction  and  we  use  this  tenninology  through  the  rest  of  this  manuscript.  TS1 
structure  is  illustrated  in  Figure  1(a).  The  H6  atom  is  in  a  bridge  structure  shifting  from  C5  to  S3. 
The  cleaving  C5-H6  bond  is  1.28A  and  the  forming  S3-H6  bond  is  1.75A.  The  cleaving  S3-C4 
bond  is  2.69A.  All  structures  are  from  B3LYP/6-31G(d,p)  determined  geometries.  The  barrier  is 
64.41  kcal/mol  is  some  20  kcal/mol  higher  than  a  conventional  6  center  retro-ene  reaction  due  to 
strain  in  the  four-center  ring  strain  of  the  transition  state. 

CH3CH2SCH2CH2C1  can  dissociate  to  CH3CH2SH  +  CH2CHC1  or  CH3CH2SC1  +  C2H4  via  a 
hydrogen  or  a  chlorine  atom  shift  (TS2B  or  TS2C),  where  transition  state  structures  are 
illustrated  in  Figures  1(c)  and  1(d),  respectively.  The  latter  reaction  (TS2C)  with  a  chlorine  shift 
has  an  ~11  kcal/mol  higher  barrier  than  TS2B  which  has  a  hydrogen  shift.  (78.61  and  67.77 
kcal/mol  barriers  for  TS2C  and  TS2B,  respectively.)  This  is  a  result  of  the  weak  S — Cl  bond 
being  fonned,  relative  to  the  S — H  bond. 

CH2C1CH2SCH2CH2C1  also  can  undergo  a  dissociation  to  CH2C1CH2SH  +  CH2CHC1  or 
CH2C1CH2SC1  +  C2H4  via  a  hydrogen  or  a  chlorine  shift  (TS3B  or  TS3C),  where  Figures  1(f) 
and  1(g)  show  transition  state  structures,  respectively.  The  optimized  geometries  of  TS3B  and 
TS3C  and  barriers  of  66.99  and  77.64  kcal/mol  in  these  reactions  are  very  similar  to  the 
optimized  geometries  of  the  respective  reactions  of  TS2B  and  TS2C  and  67.77  and  78.61 
kcal/mol  barriers  for  CH3CH2SCH2CH2C1. 

HC1  Elimination  CH3CH2SCH2CH2C1  can  dissociate  to  CH3CH2SCH=CH2  +  HC1  via  HC1 
elimination  (TS2A)  with  a  57.27  kcal/mol  barrier,  which  is  ~  10  and  21  kcal/mol  lower  than  the 
4  center  hydrogen  and  chlorine  shift  reactions  (TS2B  or  TS2C),  respectively.  The  TS2A 
structure  is  illustrated  in  Figure  1(b).  The  leaving  C4-H9  and  C5-C16  bonds  are  1.25  and  2.6lA, 
respectively. 
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CH2CICH2SCH2CH2CI  also  can  undergo  dissociation  to  CH2C1CH2SCH=CH2  +  HC1  via  HC1 
elimination  (TS3A).  The  optimized  geometry  of  TS3A  and  a  barrier  of  59.11  kcal/mol  in  this 
reaction  are  very  similar  to  the  optimized  geometry  of  TS2A  and  57.27  kcal/mol  barrier  in 
CH3CH2SCH2CH2CI  reaction  with  only  one  Cl  in  the  dithiol  molecule.  Figure  1(e)  shows 
transition  state  structure.  The  leaving  C4-H9  and  C5-C16  bonds  are  1.26  and  2.60A,  respectively. 

Bond  Cleavage  Reaction  Bond  dissociation  energies  of  CH3CH2SCH2CH3  and 
CH3CFI2SCIT2CH2CI  are  illustrated  in  Figure  5.  CFI3CH2SCFI2CFI3  can  undergo  bond  cleavage 
reactions  to  CXfiCTUS*  +  C2FI5  or  CF^CFFSCFF*  +  CFI3  with  barriers  of  73.02  or  83.95 
kcal/mol,  respectively.  CH3CH2SCH2CH2CI  can  undergo  similar  bond  cleavage  reactions  to 
CFFCFFS*  +  OH2CH2CI  or  CFFX'FFSCFF*  +  CFI2CI  with  74.42  or  83.20  kcal/mol  barriers, 
respectively.  CF^CFFSCFFCFFCl  also  can  undergo  dissociation  to  CFhCF^SCF^CF^*  +  Cl 
with  a  barrier  of  80.45  kcal/mol.  Hydrogen  BDEs  for  CH3CH2SCH2CH3  and 
CH3CH2SCH2CH2CI  were  discussed  previously. 

Thermodynamic  and  Kinetic  Analysis  of  Unimolecular  Dissociation  Reactions  Thermodynamic 
and  kinetic  analysis  on  the  unimolecular  dissociation  of  CH3CH2SCH2CH3, 
CH3CH2SCH2CH2CI,  and  CH2C1CH2SCH2CH2C1. 

Multifrequency  quantum  Rice-Ramsperger-Kassel  (QRRK)  calculations  ’  for  k(E)  and 
modified  strong  collision  analysis  of  Gilbert  et  al.  "  for  falloff  are  performed  on  this  reaction 
system  to  estimate  rate  constants  and  to  determine  important  reaction  paths  as  a  function  of 
temperature.  All  dissociation  paths  are  pressure-independent  above  0.01  atm.  Plots  of  rate 
constants  at  1  atm.  versus  1000/T  for  unimolecular  dissociations  of  CH3CH2SCH2CH3, 
CH3CH2SCH2CH2CI,  and  CH2CICH2SCH2CH2CI  are  given  in  Figure  6,  7  and  8,  respectively. 

For  CH3CH2SCH2CH3  dissociation,  bond  cleavage  reactions  to  CTfiCTXS*  +  C2H5  is  most 
important  above  800  K.  CH3CH2SH  +  C2H4  via  TS1  and  CFECFFSCFF*  +  CH3  products  have 
two  orders  of  magnitude  lower  rate  constants  than  C — S  bond  cleavage  reaction  above  1500  K. 

For  the  unimolecular  dissociations  of  CH3CH2SCH2CH2CI  and  CH2CICH2SCH2CH2CI,  HC1 
eliminations  (TS2A  or  TS3A)  are  most  important  channel  at  high  and  low  temperatures.  Carbon- 
sulfur  and  carbon-chlorine  bond  cleavage  reactions  become  important  at  higher  temperatures 
(above  1000  K),  because  of  higher  Arrhenius  pre-exponential  factors. 

Reaction  Mechanism 

Thermodynamic  data  and  kinetics  generated  in  this  work  were  combined  with  existing 

35  36  37 

mechanisms  for  combustion  chemistry  of  chlorine  ’  ,  sulfur  ,  and  nitrogen  and  hydrocarbon 
chemistry  .  The  resulting  detailed  chemical  kinetic  mechanism  models  the  complete  oxidation 
and  pyrolysis  process  for  HD  and  other  chlorosulfides. 

Sensitivity  Analysis  Using  this  detailed  kinetic  mechanism,  normalized  first  order  sensitivity 
coefficients  were  calculated  in  a  perfectly  stirred  reactor  (PSR)  code  for  a  variety  of  inlet 
compositions,  temperatures,  and  reactor  residence  times  to  determine  the  main  destruction 
pathways  for  HD. 
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Sensitivity  analysis  was  conducted  under  the  following  conditions,  chosen  to  approximate  those 
that  might  be  found  in  an  incinerator:  1)  1000  ppm  HD  in  air,  2)  a  fuel-rich  (equivalence  ratio  = 
3.0)  mixture  of  HD  and  air,  3)  1000  ppm  of  HD  in  equilibrium  fuel  lean  (equivalence  ratio  =  0.8) 
methane-air  products,  and  4)  1000  ppm  of  HD  in  fuel  lean  (equivalence  ratio  =  0.8)  methane-air 
products  with  major  species  (H2O,  O2,  CO,  CO2,  N2)  in  equilibrium  and  with  superequilibrium 
concentrations  of  H,  O,  OH,  and  Cl  radicals.  The  superequilibrium  radical  concentrations  for  the 
radicals  were  chosen  as  the  peak  values  from  a  series  of  PSR  calculations  in  which  the  residence 
time  was  varied.  The  analysis  was  conducted  for  residence  times  from  10"6  to  1000  sec.  and 
temperatures  of  1000,  1250,  1500,  1750,  and  2000  K. 

It  was  found  that  under  almost  all  conditions  that  destruction  of  the  HD  molecule  is  dominated 
by  the  HC1  elimination  pathway  described  above.  The  only  exception  to  this  is  for  HD  in  air,  at 
low  temperatures  (less  than  about  1250  K)  and  high  Cl  atom  concentrations,  where  H  abstraction 
by  Cl  is  the  dominant  pathway.  For  HD  mixed  with  methane-air  products,  HC1  elimination  was 
always  the  dominant  pathway.  Key  results  of  the  sensitivity  analysis  are  shown  in  Figures  9  and 
10. 

Summary 

Thermodynamic  properties  for  reactants,  transition  states  and  products  for  the  unimolecular 
dissociation  of  CH3CH2SCH2CH3,  CH3CH2SCH2CH2CI,  and  CH2CICH2SCH2CH2CI  are 
calculated  at  the  B3LYP/6-31G(d,p)  level.  Af//°29s  values  are  determined  with  isodesmic 
reactions.  S° 29s  and  Cp(  T)  are  determined  using  geometric  parameters  and  vibrational  frequencies 
obtained  at  the  B3LYP/6-31G(d,p).  The  hindered  rotational  contributions  to  S°29s  and  Cp(  T)  are 
calculated  by  intramolecular  torsion  potential  curves.  Hydrogen  bond  dissociation  energies  of 
CH3SH,  CH3CH2SH,  CH3SCH3,  CH3CH2SCH2CH3,  and  CH3CH2SCH2CH2CI  are  calculated  at 
the  B3LYP/6-31G(d,p).  For  CH3CH2SCH2CH3  dissociation,  bond  cleavage  reactions  to 
CFFCFFS*  +  C2H5  is  most  important  at  higher  temperature  (above  800  K).  For  the  unimolecular 
dissociations  of  CH3CH2SCH2CH2CI  and  CH2CICH2SCH2CH2CI,  HC1  eliminations  via  TS2A 
and  TS3A,  respectively,  are  the  most  important  channel  at  high  and  low  temperatures. 
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Species  Name  & 
Structure 


Bond  Length 

(A) 


Bond  Angle 
(Degree) 


Dihedral  Angle 
(Degree) 


[a]  TS1 


i-21 

r32 
i-43 
i-54 
r65 
i-75 
r85 
i-94 
rl04 
rl  12 
rl22 
rl31 
1-141 
rl51 


1.53 

1.84 

2.69 

1.40 

1.28 

1.09 

1.09 

1.08 

1.08 

1.10 

1.10 

1.10 

1.09 

1.09 


Z321 

Z432 

Z543 

Z654 

Z754 

Z854 

Z943 

Z1043 

Z1121 

Z1221 

Z1312 

Z1412 

Z1512 


114.9 

166.9 
78.1 

100.8 

117.1 

117.1 

96.6 

96.9 

111.3 

111.0 

111.0 

111.0 

111.0 


Z4321 

Z5431 

Z6543 

Z7543 

Z8543 

Z9436 

Z10436 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


96.0 

-104.3 

0.4 

-108.7 

109.2 

-121.3 

120.5 

-122.3 

120.0 

-179.0 

60.7 

-60.0 


[b]  TS2A 


r21 
r32 
i-43 
i-54 
r65 
r75 
r85 
r94 
rl04 
rl  12 
rl22 
rl31 
1-141 
rl51 


1.53 

1.85 

1.80 

1.40 

2.61 

1.08 

1.08 

1.25 

1.09 

1.09 

1.09 

1.09 

1.09 

1.09 


Z321 

Z432 

Z543 

Z654 

Z754 

Z854 

Z943 

Z1043 

Z1121 

Z1221 

Z1312 

Z1412 

Z1512 


109.4 

100.6 

119.8 
90.7 

121.1 

121.3 

109.2 
117.0 

111.3 

110.9 
109.6 

111.4 

111.5 


Z4321 

Z5432 

Z6543 

Z7543 

Z8543 

Z9432 

Z10432 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


176.0 

134.1 

-106.5 

161.8 

-16.7 

45.3 
-74.5 

-119.5 

119.4 

179.6 

60.3 
-60.8 


[c]  TS2B 


H(14) 
f  H(1 1) 

H(15Vi2,  “ 


H(1 


(9) 


r21 
r32 
i-43 
i-54 
r65 
i-75 
r85 
i-94 
rl04 
rl  12 
rl22 
rl31 
1-141 
rl51 


1.53 

1.84 

2.64 

1.40 

1.29 

1.79 

1.09 

1.08 

1.08 

1.10 

1.10 

1.10 

1.09 

1.09 


Z321 

Z432 

Z543 

Z654 

Z754 

Z854 

Z943 

Z1043 

Z1121 

Z1221 

Z1312 

Z1412 

Z1512 


114.9 

165.4 

78.4 

100.4 
118.1 
118.1 

99.9 

94.4 
111.6 
111.2 
110.7 
111.0 
111.0 


Z4321 

Z5431 

Z6543 

Z7543 

Z8543 

Z9436 

Z10436 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


95.9 

-108.9 

1.9 

-110.5 

111.1 

-122.4 

117.8 
-122.0 

119.8 
-179.4 

60.7 

-59.2 


Figure  1.  Geometries  of  Transition  States  Optimized  at  B3LYP/6-31G(d,p). 
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Species  Name  & 
_ Structure 

[d]  TS2C 


Bond  Length  Bond  Angle  Dihedral  Angle 
(A) (Degree) (Degree) 


i-21 

r32 
i-43 
i-54 
r63 
i-75 
r85 
r94 
rl04 
rl  12 
rl22 
rl31 
1-141 
rl51 


1.53 

1.84 

2.30 

1.38 

2.66 

1.08 

1.08 

1.09 

1.09 

1.10 

1.10 

1.10 

1.09 

1.09 


Z321 

Z432 

Z543 

Z632 

Z754 

Z854 

Z943 

Z1043 

Z1121 

Z1221 

Z1312 

Z1412 

Z1512 


116.2 

125.3 
95.3 

141.2 

121.4 
121.4 

97.8 

98.0 

111.1 

112.7 

110.7 
111.0 

111.2 


Z4321 

Z5431 

Z6321 

Z7543 

Z8543 

Z9432 

Z10432 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


73.2 
179.9 

-107.0 

-91.1 

91.3 
59.1 

-59.4 

-115.4 

126.8 

175.6 

55.3 
-64.8 


[e]  TS3A 


r21 

1.52 

r32 

1.85 

Z321 

108.2 

r43 

1.80 

Z432 

99.9 

r54 

1.41 

Z543 

119.4 

r65 

2.60 

Z654 

90.4 

r75 

1.08 

Z754 

120.9 

r85 

1.08 

Z854 

121.4 

r94 

1.26 

Z943 

109.3 

rl04 

1.09 

116.8 

rl  12 

1.09 

Z1043 

110.9 

rl22 

1.09 

Z1121 

110.4 

r  1 3 1 

1.81 

Z1221 

110.2 

r  1 4 1 

1.09 

Z1312 

111.8 

r  1 5 1 

1.09 

Z1412 

Z1512 

111.8 

Z4321 
Z5432 
Z6543 
Z7543 
Z8543 
Z9432 
Z 10432 
Z11213 
Z12213 
Z13123 
Z14123 
Z15123 


175.2 

134.6 
-106.3 

161.7 
-16.5 

45.4 
-74.7 

-119.6 

119.2 
-179.1 

62.5 
-60.5 


[f]  TS3B 


H(14)  H(11) 

~y.  i  _ 

H(15)  -1(12) 


Cl(7 

H(8)  L#*— 


H(10)  H(9) 


r21 

1.52 

r32 

1.84 

Z321 

111.7 

r43 

2.65 

Z432 

166.1 

r54 

1.40 

Z543 

78.4 

r65 

1.29 

Z654 

100.3 

r75 

1.78 

Z754 

118.2 

r85 

1.09 

Z854 

118.0 

i-94 

1.08 

Z943 

99.9 

rl04 

1.08 

94.1 

rl  12 

1.09 

Z1043 

111.3 

rl22 

1.09 

Z1121 

111.1 

r  1 3 1 

1.83 

Z1221 

111.1 

rl41 

1.09 

Z1312 

111.7 

r  1 5 1 

1.09 

Z1412 

Z1512 

111.7 

Z4321 

Z5431 

Z6543 

Z7543 

Z8543 

Z9436 

Z10436 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


94.0 

-112.0 

2.0 

-110.3 

111.0 

-122.6 

117.7 

-120.9 

119.9 

-179.7 

61.7 

-61.0 


Figure  1.  Geometries  of  Transition  States  Optimized  at  B3LYP/6-31G(d,p)  (continued). 
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Species  Name  & 
_ Structure 

[g]  TS3C 


CKp) 

/  i 


HjrV>  Jk.  W) 

%=e(2)T  jC£)H_(7 

JKh(i: 


Cl(15) 


12) 


ho® 

H§9) 


Bond  Length  Bond  Angle  Dihedral  Angle 
(A) (Degree) (Degree) 


r21 

1.52 

r32 

1.85 

Z321 

115.9 

i-43 

2.28 

Z432 

124.5 

i-54 

1.39 

Z543 

95.3 

r63 

2.64 

Z632 

140.5 

i-75 

1.08 

Z754 

121.5 

r85 

1.08 

Z854 

121.4 

r94 

1.09 

Z943 

98.4 

rl04 

1.09 

97.1 

rl  12 

1.10 

Z1043 

109.0 

rl22 

1.09 

Z1121 

111.9 

rl  3 1 

1.09 

Z1221 

111.5 

r  1 4 1 

1.09 

Z1312 

111.6 

rl  5 1 

1.82 

Z1412 

Z1512 

111.3 

Z4321 

Z5431 

Z6321 

Z7543 

Z8543 

Z9432 

Z10432 

Z11213 

Z12213 

Z13123 

Z14123 

Z15123 


73.1 

-179.0 

-106.5 

-91.9 

89.9 

60.3 

-58.0 

-114.9 

127.7 

173.1 

49.8 

-68.9 


Figure  1.  Geometries  of  Transition  States  Optimized  at  B3LYP/6-31G(d,p)  (continued). 
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Rotational  Barrier  (kcal/mol)  Rotational  Barrier  (kcal/mol) 


0  60  120  180  240  300  360  0  60  120  180  240  300  360 

HC--CS  Torsion  Angle  +  60.0  (degree)  CC--SC  Torsion  Angle  -  177.2  (degree) 


CS-CC  Torsion  Angle  -  80.2  (degree)  SC--CC1  Torsional  Angle  -  180.0  (degree) 


Figure  2.  Potential  barriers  for  4  internal  rotations  about  CH3CH2SCH2CH2CI.  Points  are 
calculated  values  at  the  B3LYP/6-31G(d,p)  level.  Lines  are  Fourier  expansion. 
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Figure  3  Potential  barriers  for  2  internal  rotations  about  CH2CICH2SCH2CH2CI.  Points  are 
calculated  values  at  the  B3LYP/6-31G(d,p)  level.  Lines  are  Fourier  expansions. 
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1)CH3CH2SCH2CH3 


CH3CH2SCH2C.H2  +  H  (BDE=98.22) 


CH3CH2SC.HCH3  +  H  (BDE=89.90) 

1 

CH3CH2SC.H2  +  CH3  (BDE=83.95) 

TS1  Ea=64.41 
(H  shift) 

CH3CH2SCH2CH3  /""\ 

1 

CH3CH2S.  +  C2H5  (BDE=73.02) 

CH3CH2SH  +  CH2=CH2  (0.48) 

(-20.08)  I 

. $---  dH°(rxn)=20.56 

2)  CH3CH2SCH2CH2CI 


CH3CH2SCH2C.HCI  +  H  (BDE=96.98) 


TS2C  Ea=78.61 
(Cl  shift) 


TS2B  Ea=67.77 
(H  shift) 

TS2A  Ea=57.27 
(HCI  Elim.) 

CH3CH2SCH2CH2CI 
(-26.45) 


CH3CH2SC.HCH2CI  +  H  (BDE=91.31) 


CH3CH2SC.H2  +  C.H2CI  (BDE=83.20) 


CH3CH2SCH2C.H2  +  Cl  (BDE=80.45) 


CH3CH2S.  +  C.H2CH2CI  (BDE=74.42) 


CH3CH2SCI  +  CH2=CH2  (-0.41) 


CH3CH2SH  +  CH2=CHCI  (-7.55) 


dH°(rxn)=26.04 

dH°(rxn)=18.90 


CH3CH2SCH=CH2  +  HCI  (-11.47) 
. -t-  dH°(rxn)=14.98 


3)  CH2CICH2SCH2CH2CI 

TS3C  Ea=77.64 
(Cl  shift) 

TS3B  Ea=66.99  * 

(H  shift) 

TS3A  Ea=59.11t  ii, 

(HCI  Elim.) 

CH2CICH2SCH2CH2CI 
(-31.12) 


CH2CICH2S.  +  C.H2CH2CI  (BDE=74.12) 


CH2CICH2SCI  +  CH2=CH2  (-4.21) 


1 


CH2CICH2SH  +  CH2=CHCI  (-12.46) 
CH2CICH2SCH=CH2  +  HCI  (-16.43) 


\XdH°(rxn)=18.66 
X  dH°(rxn)=26.91 


dH°(rxn)=14.69 


Figure  4.  Potential  Energy  Diagram  (Units:  kcal/mol). 
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(a)  CH3CH2SCH2CH3 


(b)  CH3CH2SCH2CH2CI 


Figure  5.  Bond  dissociation  energies3  (kcal/mol)  of  CH3CH2SCH2CH3b  and 
CH3CH2SCH2CH2C1c. 

a  Bond  energies  are  calculated  in  this  work.  (Table  5) 

b  CH3CH2SCH2 — CH3:  AfH° 298  of  [C2H5SCH3]  =  -14.24  ±  0.26  (ref.  7),  [CH3]  =  34.82  ±  0.2 
(ref.  8)  and  BDE  of  [CH3SCH2 — -H]  =  95.39  (in  this  study)  are  used. 
c  CH3CH2S— CH2CH2CI:  Af//° 29g  of  [C*H2CH2C1]  =  23.83  (ref.  36)  is  used. 
c  CH3CH2SCH2 — CH2CI:  Af//°298  of  [OH2Cl]  =  27.7  ±  2.0  (ref.  37)  is  used. 
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Figure  6.  CH3CH2SCH2CH3  Dissociation  constant  vs.  1000/T  at  1  atm. 


Figure  7.  CH3CH2SCH2CH2CI  Dissociation  constant  vs.  1000/T  at  1  atm. 
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8 


P  =  1  atm. 


Figure  8.  CFI^CICFFSCFFCFFCI  Dissociation  constant  vs.  1000/T  at  1  atm. 


Fraction  of  HD  Remaining 
residence  time  decreasing . > 


Fraction  of  HD  Remaining 
residence  time  decreasing . > 

Figure  9.  Calculated  first  order  sensitivity  coefficients  for  FID  destruction  in  a  perfectly 
stirred  reactor  for  1000  ppm  of  HD  in  air,  P  =  1  atm. 
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1.E+01 


1.E+00 

>>  1.E-01 
> 

|  1.E-02 

c 

0) 

1.E-03 

1.E-04 


1.E-05 


js  ♦  >:  ♦  >:  ♦  1  ♦  kui  ♦  ■>: 


(a)  HCI  Elimination 


- 1000  K 
- 1250  K 
1500  K 
1750  K 
2000  K 


1.E-12  1.E-10  1.E-08  1.E-06  1.E-04  1.E-02  1.E+00 

Fraction  of  HD  Remaining 
residence  time  decreasing . > 


Figure  10.  Calculated  first  order  sensitivity  coefficients  for  FID  destruction  in  a  perfectly 
stirred  reactor  for  1000  ppm  of  HD  in  lean  (equivalence  ratio  =  0.8)  CH4-air  products  with 
superequilibrium  radicals,  P  =  1  atm. 
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TABLE  1:  Comparison  of  Enthalpies  of  Formation  at  B3LYP/6-31G(d,p)  with  Literature 
Values  (Units  in  kcal/mol). 


Working  Reactions 


B3LYP/6- 

31G(d,p) 


A//f°  298 


\xnfl  298 

Af H°  298 

Avg. 

(B3LYP)a 

Literature 

CH3SH  +  CH4  <=>  CH3CH3  +  h2s 

-1.26 

-5.99 

-6.22  ±  0.68 

-5.47  ±0.14(7] 

CH3SH  +  CH3OH  <=>  ch3ch2oh  +  h2s 

-6.47 

-6.44 

-11.07  ± 

CH3CH2SH  +  CH4  <=>  CH3SH  +  CH3CH3 

2.52 

-10.34 

-10.82  ±2.93 

0.14(7] 

CH3CH2SH  +  CH3OH  <=>  CH,SH  +  CH3CH2OH 

-2.69 

-10.79 

ch3ch2sh  +  CH3CH2OH  <=>  CH3SH  + 
(CH3)2CHOH 

-3.17 

-11.32 

CH3SCH3  +  CH4  <=>  CH3SH  +  CH3CH3 

0.87 

-8.69 

-9.17  ±  2.93 

-8.96  ±0.12(7] 

CH3SCH3  +  CH3OH  <=>  CH3SH  +  CH3CH2OH 

-4.33 

-9.15 

CH3SCH3  +  CH3CH2OH  <=>  CH3SH  + 

(CH3)2CHOH 

-4.82 

-9.67 

36.36  ± 

C.H2SH  +  CH3OH  <=>  CH,SH  +  C.H2OH 

-1.13 

39.85 

38.72  ±2.12 

2.00(14] 

c.h2sh  +  ch4  <=>  c.h2ch,  +  h2s 

3.43 

38.36 

C.HjSH  +  CH3OH  <=>  C.H2CH2OH  +  H2S 

-0.38 

37.94 

ch3s.  +  ch,ch2oh<=>  ch3sh  +  ch3ch2o. 

15.97 

30.83 

31.39  ±2.29 

30.83  ± 
1.90(18] 

CH3S.  +  CH3OH  <=>  CH3SH  +  ch3o. 

14.85 

31.94 

ch3ch2s.  +  ch3ch2oh  <=>  ch3ch2sh  + 
ch3ch2o. 

15.73 

25.47 

25.93  ±  2.09 

24.14  ±0.96(3] 

CH3CH2S.  +  CH3OH  <=>  CHjSH  +  CH3CH20. 

13.05 

25.74 

CH3CH2S.  +  CH3OH  <=>  CH3CH2SH  +  CH3O. 

14.62 

26.57 

c.h2ch2sh  +  CH3SH  <=>  ch3ch2sh  +  c.h2sh 

-3.57 

34.33 

35.60  ±4.06 

C.H2CH2SH  +  CH3CH3  <=>  ch3ch2sh  +  c.h2ch3 

1.11 

36.86 

CHjC.HSH  +  CH3SH  <=>  CH3CH2SH  +  C.H2SH 

2.67 

28.09 

29.35  ±  4.06 

CHjC.HSH  +  CH3CH3  <=>  CH3CH2SH  +  C.H2CH3 

7.36 

30.61 

C.H2SCH3  +  CH3OCH3  <=>  CH3SCH3  +  CH3OC.H2 

0.68 

34.45 

34.33  ±  1.23 

32.27  ± 
0.72(20,21] 

C.H2SCH3  +  CH3CH2OH  <=>  CH3SCH3  + 
C.H2CH2OH 

7.30 

34.21 

a  Uncertainties  are  sum  of  the  standard  deviation  for  B3LYP/6-31G(d,p)  level  of  theory  and  maximum 
cumulative  uncertainties  from  reference  species. 


21 


TABLE  2:  Enthalpies  of  Formation  for  Reference  Molecules  in  the  Isodesmic  Reactions 


Compounds 

AfH°  29g  (kcal/mol) 

Compounds 

Af77°  29g  (kcal/mol) 

CH3SH 

-5.47  ±0.14  [7] 

ch4 

-17.89  ±  0.07a  [8] 

CH3CH3 

-20.24  ±0.10a  [9] 

h2s 

-4.90  ±0.19  [10] 

CH3CH2SH 

-1 1.07  ±  0.14  [7] 

CH3OH 

-48.16  ±0.07  [11] 

CH3CH2OH 

-56.17  ±0.10a  [12] 

(CH3)2CHOH 

-65.19  ±2.2  [13] 

CH3SCH3 

-8.96  ±0.12  [7] 

c.h2sh 

36.36  ±2.0  [14] 

c.h2oh 

-3.97  ±0.22  [15] 

C.H2CH3 

28.80  ±0.50  [16] 

C.H2CH20H 

-5.70  ±0.85  [17] 

CH3S. 

30.83  ±  1.90  [18] 

CH3CH20. 

-3.90  ±  1.27  [15] 

CH3O. 

4.10  ±  1.0  [19] 

ch3ch2s. 

24.14  ±0.96  [3] 

C.H2SCH3 

32.27  ±0.72  [20,21] 

ch3ch2sch2ch3 

-20.08  ±  0.24  [7] 

CH3CH0CI 

-26.80  ±  0.26a  [8] 

CH0C1CH2C1 

-30.33  ±0.67  [11] 

CH3OCH3 

-43.99  ±0.12  [11] 

CH30C.H2 

0.1  [22] 

CH3C.HOH 

-13.34  ±0.84  [15] 

a  The  uncertainties  are  evaluated  from  reference  [11]. 
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TABLE  3:  Reaction  Enthalpies  and  Enthalpies  of  Formation  in  the  Isodesmic  Reactions 
(Units  in  kcal/mol). 


Working  Reactions 

Arxn H0 298 

Af H°  298 

B3LYP 

KMLYP 

B3LYP 

KMLYP 

C.H2CH2SH  +  CH3SH  <=>  CH3CH2SH  +  C.H2SH 

-3.57 

-4.45 

34.33 

35.21 

CH3C.HSH  +  CH3SH  <=>  CH3CH2SH  +  C.H2SH 

2.67 

1.82 

28.09 

28.94 

ch3ch2sch2c.h2  +  ch3ch2sh  <=>  ch3ch2sch2ch3  +  c.h2ch2sh 

-0.72 

-0.49 

26.04 

26.69 

CH3CH2SC.HCH3  +  CH3CH2SH  <=>  CH3CH2SCH2CH3  +  CHjC.HSH 

1.36 

1.21 

17.72 

18.72 

CH3CH2SCH2CH2C1  +  ch3ch3  <=>  CH3CH2SCH2CH3  +  CH3CH2C1 

-0.19 

-0.23 

-26.45 

-26.41 

ch2cich2sch2ch2ci  +  ch3ch3  <=>  CH3CH2SCH2CH3  +  CH2C1CH2C1 

0.95 

1.18 

-31.12 

-31.35 

ch3ch2sch2c.hci  +  CH3CH2OH  <=>  CH3CH2SCH2CH2C1  + 
C.H2CH2OH 

5.71 

4.70 

18.32 

19.36 

ch3ch2sch2c.hci  +  CH3OCH3  <=>  CH3CH2SCH2CH2C1  +  CH3OC.H2 

-0.91 

-1.69 

18.55 

19.38 

Average2  for  CH3CH2SCH2C.HC1 

CH3CH2SC.HCH2C1  +  CH3CH2OH  <=>  CH3CH2SCH2CH2C1  + 

CHjC.HOH 

2.80 

0.99 

13.58 

18.43  ± 

1.71 

(B3LYP) 
19.37  ± 

1.56 

(KMLYP) 

15.43 

CH3CH2SC.HCH2C1  +  CH3OCH3  <=>  CH3CH2SCH2CH2C1  +  CH3OC.H2 

5.70 

3.06 

11.95 

14.62 

Average2  for  CH3CH2SC.HCH2Clb 

12.76  ± 
2.69 
(B3LYP) 
15.03  ± 
2.11 
(KMLYP) 

2  Uncertainties  are  sum  of  the  standard  deviation  for  each  calculation  level  of  theory  and  maximum 
cumulative  uncertainties  from  reference  species. 

b  The  values  between  B3LYP  and  KMLYP  have  2.3  kcal/mol  difference,  so  further  calculations  were 
performed.  12.71  ±  2.18  (B3LYP/6-311G(d,p))  and  12.85  ±  1.65  (CBSQ),  so  B3LYP  value  is 
recommended. 
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TABLE  4:  Comparison  of  Activation  Energies  (Units  in  kcal/mol). 


B3LYP/6- 

31G(d,p) 

KMLYP/6- 
31  lG(d,p) 

CBS-Q 

CBS- 

QB3 

<H  Shift> 

TS1 

61.51 

73.42 

65.46 

64.41 

TS2B 

65.62 

78.08 

68.62 

67.77 

TS3B 

64.15 

76.93 

66.81 

66.99 

<C1  Shift> 

TS2C 

73.54 

95.29 

80.89 

78.61 

TS3C 

73.31 

94.80 

79.57 

77.64 

<HC1  Elimination> 

TS2A 

52.25 

59.04 

58.06 

57.27 

TS3A 

54.27 

61.42 

60.08 

59.11 
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TABLE  5:  Ideal  Gas  Phase  Thermodynamic  Properties  Obtained  by  B3LYP/6-31G(d,p) 
Calculation3 


Species 
(  s,  e,  OI  )8 

Af H° 298b 

rtO  C 

^  298 

Q?3ooC 

Q?400 

Q?500 

Q?600 

Cp  800 

Q?1000 

Q?1500 

REACTANTS 

CH3CH2SCH2CH3 

TVRd 

66.19 

19.36 

25.96 

32.18 

37.59 

46.23 

52.73 

62.98 

(C-CSCC) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCC) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CC) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCSC-C) 

I.  R. 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(18,0,1) 

Total 

-20.08 

89.97 

27.02 

33.34 

39.18 

44.17 

52.11 

58.13 

67.70 

CH3CH2SCH2CH2C1 

TVRd 

75.47 

22.00 

28.87 

35.10 

40.38 

48.62 

54.70 

64.16 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CCC1) 

I.  R. 

7.44 

2.27 

1.92 

1.73 

1.61 

1.47 

1.37 

1.23 

(CCSC-CC1) 

I.  R. 

5.73 

3.67 

3.77 

3.49 

3.12 

2.51 

2.11 

1.59 

(3,0,1) 

Total 

-26.45 

100.53 

31.77 

38.25 

43.82 

48.40 

55.54 

60.88 

69.34 

CH2C1CH2SCH2CH2C1 

TVRd 

82.45 

24.60 

31.74 

38.00 

43.16 

51.01 

56.66 

65.34 

(C1C-CSCCC1) 

I.  R.f 

5.95 

3.75 

3.71 

3.38 

3.01 

2.43 

2.06 

1.57 

(C1CC-SCCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(C1CCS-CCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(C1CCSC-CC1) 

I.  R. 

5.95 

3.75 

3.71 

3.38 

3.01 

2.43 

2.06 

1.57 

(1,0,1) 

Total 

-31.12 

110.07 

36.38 

42.70 

47.98 

52.22 

58.71 

63.48 

70.92 

TRANSITION  STATES 

TS1 

TVRd 

74.88 

24.34 

31.13 

37.25 

42.45 

50.59 

56.62 

66.06 

(C-CSCC) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCC) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(3,0,1) 

Total 

44.33 

86.77 

28.17 

34.82 

40.75 

45.74 

53.53 

59.32 

68.42 

TS2A 

TVRd 

78.71 

25.12 

31.93 

37.96 

43.01 

50.79 

56.45 

65.20 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CCC1) 

I.  R. 

7.44 

2.27 

1.92 

1.73 

1.61 

1.47 

1.37 

1.23 

(3,0,1) 

Total 

30.82 

98.04 

31.22 

37.54 

43.19 

47.91 

55.20 

60.52 

68.79 

TS2B 

TVRd 

81.67 

27.31 

34.28 

40.34 

45.37 

53.08 

58.68 

67.33 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(3,0,1) 

Total 

41.32 

93.56 

31.14 

37.97 

43.84 

48.66 

56.02 

61.38 

69.69 

TS2C 

TVRd 

85.11 

27.96 

34.32 

39.94 

44.67 

52.09 

57.66 

66.55 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(3,0,1) 

Total 

52.16 

97.00 

31.79 

38.01 

43.44 

47.96 

55.03 

60.36 

68.91 

TS3A 

TVRd 

87.38 

27.74 

34.77 

40.82 

45.77 

53.16 

58.42 

66.40 

(C1C-CSCCC1) 

I.  R.f 

5.95 

3.75 

3.71 

3.38 

3.01 

2.43 

2.06 

1.57 

(C1CC-SCCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(C1CCS-CCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(1,0,1) 

Total 

27.99 

109.05 

35.77 

42.02 

47.42 

51.82 

58.43 

63.18 

70.41 

TS3B 

TVRd 

90.45 

29.93 

37.10 

43.17 

48.09 

55.42 

60.61 

68.49 

(C1C-CSCCC1) 

I.  R.f 

5.95 

3.75 

3.71 

3.38 

3.01 

2.43 

2.06 

1.57 

(C1CC-SCCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(1,0,1) 

Total 

35.87 

104.26 

35.82 

42.58 

48.16 

52.62 

59.27 

64.02 

71.28 

TS3C 

TVRd 

92.58 

30.43 

37.08 

42.74 

47.38 

54.44 

59.61 

67.74 

(C1C-CSCCC1) 

I.  R.f 

5.95 

3.75 

3.71 

3.38 

3.01 

2.43 

2.06 

1.57 
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TABLE  5:  (Continued)  Ideal  Gas  Phase  Thermodynamic  Properties  Obtained  by  B3LYP/6-31G(d,p)  Calculation3 


(C1CC-SCCC1) 

I.  R. 

7.86 

2.14 

1.77 

1.61 

1.52 

1.42 

1.35 

1.22 

(1,0,1) 

Total 

46.52 

106.39 

36.32 

42.56 

47.73 

51.91 

58.29 

63.02 

70.53 

PRODUCTS 

CH3CH2SCH2C.H2 

TVRd 

71.66 

20.15 

26.31 

31.88 

36.61 

44.09 

49.72 

58.66 

(C-CSCC) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCC) 

I.R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CC) 

I.R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCSC-C) 

I.R. 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(3, 1/2,1) 

Total 

26.04 

95.44 

27.81 

33.69 

38.88 

43.19 

49.97 

55.12 

63.38 

ch3ch2sc.hch3 

TVRd 

69.19 

19.66 

25.69 

31.27 

36.08 

43.74 

49.50 

58.60 

(C-CSCC) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCC) 

I.R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CC) 

I.R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCSC-C) 

I.R. 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(9, 1/2,1) 

Total 

17.72 

92.97 

27.32 

33.07 

38.27 

42.66 

49.62 

54.90 

63.32 

CH3CH2SCH2C.HC1 

TVRd 

77.13 

22.39 

28.71 

34.29 

38.96 

46.20 

51.53 

59.83 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CCC1) 

I.  R. 

7.44 

2.27 

1.92 

1.73 

1.61 

1.47 

1.37 

1.23 

(CCSC-CC1) 

I.R. 

5.73 

3.67 

3.77 

3.49 

3.12 

2.51 

2.11 

1.59 

(3,1/2, 1) 

Total 

18.43 

102.19 

32.16 

38.09 

43.01 

46.98 

53.12 

57.71 

65.01 

CH3CH2SC.HCH2C1 

TVRd 

76.51 

22.33 

28.62 

34.18 

38.84 

46.07 

51.39 

59.71 

(C-CSCCC1) 

I.  R.1 

4.40 

2.11 

2.12 

2.02 

1.88 

1.63 

1.46 

1.23 

(CC-SCCC1) 

I.  R. 

7.49 

1.72 

1.57 

1.48 

1.41 

1.31 

1.24 

1.13 

(CCS-CCC1) 

I.  R. 

7.44 

2.27 

1.92 

1.73 

1.61 

1.47 

1.37 

1.23 

(CCSC-CC1) 

I.R. 

5.73 

3.67 

3.77 

3.49 

3.12 

2.51 

2.11 

1.59 

(3,1/2, 1) 

Total 

12.76 

101.57 

32.10 

38.00 

42.90 

46.86 

52.99 

57.57 

64.89 

ch3ch2sh 

TVRd 

62.12 

13.43 

17.03 

20.40 

23.34 

28.11 

31.74 

37.51 

(C-CSH) 

I.  R.1 

4.27 

2.10 

2.16 

2.09 

1.96 

1.71 

1.53 

1.27 

(CC-SH) 

I.R. 

4.57 

1.94 

1.64 

1.44 

1.32 

1.19 

1.12 

1.05 

(3,0,1) 

Total 

-10.82 

70.96 

17.47 

20.83 

23.93 

26.62 

31.01 

34.39 

39.83 

ch3ch2s. 

TVRd 

63.58 

13.65 

16.85 

19.82 

22.39 

26.53 

29.67 

34.64 

(C-CS) 

I.  R.1 

4.50 

2.09 

2.06 

1.94 

1.79 

1.55 

1.40 

1.20 

(3, 1/2,1) 

Total 

25.93 

68.08 

15.74 

18.91 

21.76 

24.18 

28.08 

31.07 

35.84 

CH3CH2SC1 

TVRd 

68.19 

15.91 

19.55 

22.82 

25.61 

30.00 

33.28 

38.42 

(C-CSCl) 

I.  R.1 

4.40 

2.11 

2.13 

2.03 

1.89 

1.64 

1.46 

1.23 

(CC-SC1) 

I.R. 

6.69 

2.21 

2.07 

1.92 

1.78 

1.56 

1.42 

1.22 

(3,0,1) 

Total 

-12.95 

79.28 

20.23 

23.75 

26.77 

29.28 

33.20 

36.16 

40.87 

ch3ch2schch2 

TVRd 

69.15 

19.38 

25.37 

30.70 

35.20 

42.23 

47.47 

55.73 

(C-CSC=C) 

I.  R.1 

4.48 

2.14 

2.11 

1.97 

1.82 

1.57 

1.40 

1.20 

(CC-SC=C) 

I.R. 

6.85 

2.13 

1.90 

1.75 

1.65 

1.49 

1.38 

1.22 

(CCS-C=C) 

I.R. 

6.17 

2.24 

2.09 

1.93 

1.78 

1.55 

1.40 

1.20 

(3,0,1) 

Total 

10.59 

86.65 

25.89 

31.47 

36.35 

40.45 

46.84 

51.65 

59.35 

CH2C1CH2SC1 

TVRd 

76.25 

18.46 

22.34 

25.63 

28.31 

32.34 

35.21 

39.60 

(C1C-CSC1) 

I.  R.1 

5.72 

3.26 

3.38 

3.30 

3.10 

2.65 

2.27 

1.72 

(C1CC-SC1) 

I.R. 

6.89 

2.39 

2.43 

2.31 

2.14 

1.83 

1.61 

1.32 

(1,0,1) 

Total 

-16.75 

88.86 

24.11 

28.15 

31.24 

33.55 

36.82 

39.09 

42.64 

a  :  Thermodynamic  properties  are  referred  to  a  standard  state  of  an  ideal  gas  of  pure  enantiometer  at  1  atm. 
b  :  Units  in  kcal/mol  c  :  Units  in  cal/mol-K 

d  :  Sum  of  contributions  from  translations,  vibrations,  and  external  rotations, 
f :  I.  R.  represents  contribution  from  internal  rotation 

g  :  Symmetry  number,  optical  isomer  and  electronic  spin  are  taken  into  account,  -Rln(s),  Rln2,  Rln2,  respectively. 

s  =  number  of  symmetry,  e  =  electronic  spin,  OI  =  number  of  optical  isomer 
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TABLE  6:  Bond  energies  at  298K  (kcal/mol). 


Reaction  series 

Bond  energies 

B3LYP3 

KMLYPb 

Literaturec 

CH3SH  =>  CH3S.  +  H. 

88.96  ±2.43 

- 

88.40  ±2.04[7, 18] 

CH3SH  =>  C.H2SH  +  H. 

96.29  ±2.26 

- 

93.93  ±2.14[7, 14] 

CH3CH2SH  =>  CH3CH2S.  +  H. 

89.10  ±2.23 

- 

87.31  ±  1.10[7,3] 

CH3CH2SH  =>  C.H2CH2SH  +  H. 

98.77  ±4.20 

98.38  ±4.20 

n/a 

CH3CH2SH  =>  CH3C.HSH  +  H. 

92.52  ±4.20 

92.11  ±4.20 

n/a 

CH3SCH3  =>  C.H2SCH3  +  H. 

95.39  ±  1.35 

- 

93.33  ±0.84[7,20,21] 

ch3ch2sch2ch3  =>ch3ch2sch2c.h2  +  H. 

98.22  ±2.90 

98.87  ±2.90 

n/a 

ch3ch2sch2ch3  =>ch3ch2sc.hch3  +  H. 

89.90  ±2.90 

90.90  ±2.90 

n/a 

CH3CH2SCH2CH2C1  =>  CH3CH2SCH2C.HC1  +  H. 

96.98  ±  1.71 

97.88  ±  1.56 

n/a 

CH3CH2SCH2CH2C1  =>  CH3CH2SC.HCH2C1  +  H. 

91.31  ±2.69 

93.53  ±2.11 

n/a 

a  Uncertainties  are  sum  of  the  standard  deviation  for  B3LYP/6-3  lG(d,p)  level  of  theory  and  maximum 
cumulative  uncertainties  from  reference  species. 

b  Uncertainties  are  sum  of  the  standard  deviation  for  KMLYP/6-3 1  lG(d,p)  level  of  theory  and  maximum 
cumulative  uncertainties  from  reference  species. 

c  Uncertainties  are  sum  of  maximum  cumulative  uncertainties  from  reference  species. 
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